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spectrum of cosmic rays 

J.C. Arteaga-Velazquez/'Q W.D. Apel,^ K. Bekk,^ M. Bertaina,^ J. Blumer,^'! H. Bozdog,^ I.M. Brancus," 
P. Buchholz,^ E. Cantoni,^''^ A. Chiavassa,^ F. Cossavella/'[l| K. Daumiller,^ V. de Souza,i'0F. Di Pierro,^ 
P. Doll,2 R. Engel,2 J. Engler,^ M. Finger,^ D. Fuhrmann/ P.L. Ghia,^ H.J. Gils,^ R. Glasstetter/ 
C. Grupen,^ A. Haungs,^ D. Heck,^ J.R. H6randel/'| T. Huege,^ P.G. Isar,^ K.-H. Kampert/ D. Kang,i 
D. Kickelbick,^ H.O. Klages,^ K. Link/ P. Luczak,® M. Ludwig,i H.J. Mathes,^ H.J. Mayer,^ M. Melissas/ 

J. Milke,2 B. Mitrica,^ C. Morello,^ G. Navarra,^.^ S. Nehls,^ J. Oehlschlager,^ S. Ostapchenko,^^ 
S. Over,5 N. Palmieri,i M. Petcu,^ T. Pierog,^ H. Rebel,^ M. Roth,^ H. Schieler,^ F. Schroder,^ O. Sinia,^ 
G. Toma,"* G.C. Trinchero,*^ H. Ulrich,^ A. Weindl,^ J. Wochele,^ M. Wommer,^ and J. Zabierowski*^ 

'Institut filr Experimentelle Kernphysik, Karlsruher Institut fiir Technologic - Campus Slid, Germany 
'^Institut fiir Kernphysik, Karlsruher Institut fiir Technologic - Campus Nord, Cermany 
Dipartimento di Fisica Generale dell' Universita Torino, Italy 
'''National Institute of Physics and Nuclear Engineering, Bucharest, Romania 
^' Fachbereich Physik, Universitat Siegen, Germany 
^Istituto di Fisica dello Spazio Interplanetario , INAF Torino, Italy 
' Fachbereich Physik, Universitat Wuppertal, Germany 
^ Soltan Institute for Nuclear Studies, lodz, Poland 
'^Department of Physics, University of Bucharest, Bucharest, Romania 

The all-particle energy spectrum as measured by the KASCADE-Grande experiment for E — 
10^^ - 10^* eV is presented within the framework of the QGSJET II/FLUKA hadronic interaction 
models. Three different methods were applied based on the muon size and the total number of 
charged particles individually and in combination. From the study it is found that the spectrum 
cannot be completely described by a smooth power law due to the presence of characteristic features. 



INTRODUCTION 



Powerful cosmic ray accelerators are hidden inside 
our own galaxy and in the heart of extragalactic ob- 
jects, but their precise location and exact working 
mechanism are still subject of intense debate. Cos- 
mic ray research provides important clues through the 
measurement of the energy spectrum, the arrival di- 
rection and the composition of such enigmatic parti- 
cles by means of the extensive air showers (EAS) that 
cosmic rays produce in the Earth atmosphere. Along 
the past years, dedicated experiments have shown that 
the energy spectrum of cosmic rays extends from a 
few MeV up to 10^" eV, has a striking power-law be- 
havior with a spectral index 7 w —2.7 and exhibits 
some intriguing structures around 3 — 5-10^^ and 
4 — 10 ■ 10^^ eV, better known as the knee and the 
ankle, where the spectral index of the particle flux de- 
creases and grows, respectively pj. The precise inter- 



pretation of these features are not yet clear, but they 
could arise as a consequence of an interplay among 
different factors such as the loss of efhciency of the 
galactic accelerators, the dominance of a new compo- 
nent of extragalactic origin at the high-energy regime 
and the propagation effects of cosmic rays through 
space 2, 3]. To test different scenarios, accurate and 
high-statistics measurements are needed between 10^^ 
and 10^^ eV. Although observations of this kind have 
been already performed around the knee, the picture 
is not yet complete due to the lack of quality EAS 
data at higher energies: E = 10^^ — 10^^ eV. The goal 
of the KASCADE-Grande experiment is to close this 
gap between the knee and the ultrahigh-energy region 
with accurate EAS measurements. 
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KASCADE-Grande is a ground-based EAS array 
integrated by different detection systems to measure 
and separate the muon and electromagnetic compo- 
nents of the EAS Once combined, both observ- 
ables become a powerful tool to reconstruct the energy 
spectrum and to estimate the composition of primary 
cosmic ray events. In the following lines, the energy 
spectrum of cosmic rays for E = 10^^ — 10^^ eV as mea- 
sured by KASCADE-Grande will be presented and the 
reconstruction method will be described. First, a brief 
introduction to the characteristics and performance of 
the experiment will be given. 
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TABLE I: Some components of KASCADE-Grande: 
total sensitive areas and threshold kinematic energies for 
vertically incident particles are presented. MTD refers to 
the muon tracking detectors 



Detector 



Particle Area(m'^) Threshold 



Grande array 
(plastic scintillators) 
Piccolo array 

(plastic scintillators) Charged 
KASCADE array 
(liquid scintillators) 
KASCADE array 
(shielded plast. scint.) 
MTD (streamer tubes) 



Charged 370 



e/7 



80 

490 

622 
4 X 128 



3 MeV 

3 MeV 

5 MeV 

230 MeV 
800 MeV 



II. THE KASCADE-GRANDE EXPERIMENT 
A. Description 

The KASCADE-grande detector (located at 49.1° 
N, 8.4° E, 110 m a.s.l.) is the successor of the KAS- 
CADE experiment and incorporates the original 
electromagnetic detectors and muon devices of KAS- 
CADE to a bigger system of detectors, called Grande, 
composed of a 700 x 700 m^ array with 37 x 10 m^ 
scintillator stations regularly spaced by an average 
distance of 137 m ^J. To coordinate the trigger of 
KASCADE with EAS events where the core is found 
inside Grande a smaller array, named Piccolo, was 
added. The main characteristics of the KASCADE- 
Grande detectors are displayed in table U and the lay- 
out of the experiment is presented in figure [T] 

The Grande array is used to sample the density of 
charged particles of the shower front at ground level 
and to measure the particle arrival times of the EAS. 
The core position, the number of charged particles 
(Nch) and the arrival direction of the shower are ex- 
tracted from the Grande data through an iterative 
fit and a careful modeling of the EAS front. To re- 
construct the arrival direction a fit is applied as- 
suming a curved shower front as suggested by COR- 
SIKA/QGSJET II simulations. On the other hand, to 
obtain the core position and the shower size a max- 
imum likelihood procedure was employed by using a 
modified NKG lateral distribution function as a refer- 
ence model \4^. 

An important component of this experiment is the 
KASCADE muon array, composed by 192 x 3.2 m^ 
shielded scintillator detectors, which are sensitive to 
muons with threshold energy above 230 MeV for ver- 
tical incidence. With the Grande information and the 
measurements of the shielded array of the lateral dis- 
tribution of muons in the shower front, the muon size 
{Nf^) is reconstructed event-by-event at KASCADE- 



Grande. The procedure involves the maximum log- 
likelihood technique along with a Lagutin-Raikin dis- 
tribution function (^]. 



B. Accuracy of EAS Reconstruction 

Systematic uncertainties for the core, Nch and ar- 
rival direction of the EAS were studied directly Q by 
comparing the results of the Grande and KASCADE 
reconstructions, which work independently, for a sub- 
set of data with cores located inside a common area 
and shower sizes in the interval logj^Q Nch = 5.8 — 7.2. 
EAS core positions and arrival directions were found 
to be estimated by Grande with an accuracy of the or- 
der of 6 m and 0.8°, respectively. The same evaluation 
showed that the systematic uncertainties of Nch at 
Grande are < 5%. Those values are in full agreement 
with expectations from Monte Carlo simulations. 

The same kind of study cannot be applied to the 
muon size since shielded detectors are only available 
for KASCADE. In this case, Monte Carlo simula- 
tions have to be used. The results show that, for 
Nch > lO**'^, the Nf^ systematic uncertainties are 
< 25%, decreasing for high energies (see figure [U left), 
and exhibit a dependence with the core position and 
the zenith angle (see, for example, figure [5J right). 
Since the behavior of the muon uncertainties are well 
understood, a muon correction function can be con- 
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FIG. 1: The KASCADE-Grande experiment. Small 
squares represent the Grande stations. The Kascade ar- 
ray is seen at the upper right hand of the figure. KAS- 
CADE detectors are arranged in 16 clusters. The outer 
12 clusters contain the muon detectors. The dotted line 
encloses the fiducial area selected for the present analy- 
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FIG. 2: Muon systematic uncertainties vs A'^ (left) and the distance to the KASCADE center (right). The situations 
before (soUd points) and after (open circles) applying the muon correction fimction are compared. Vertical lines and 
grey bands represent the one-sigma errors, respectively. Plots were calculated with CORSIKA/QGSJET II for a mixed 
composition. Dotted lines correspond to pure proton and iron nuclei after using the muon correction function (upper 
and lower dashed lines, respectively). 



structed and applied to the muon data. This function 
was built from Monte Carlo simulations for a mixed 
composition assumption and was parameterized as a 
function of the muon size, the EAS arrival direction 
and the distance to the KASCADE center. The sys- 
tematic uncertainties of the corrected muon number 
are < 15% at threshold and quickly decrease for higher 
muon numbers as seen in figure [21 



III. THE ENERGY SPECTRUM 

A. Data sample 

In order to reduce the influence of systematic uncer- 
tainties in the data, a fiducial area located at the cen- 
ter of KASCADE-Grande (as shown in figure [T]) was 
selected and only events with zenith angles < 40°, 
which passed the KASCADE-Grande reconstruction 
without failures were considered. Additionally, sev- 
eral experimental cuts were imposed, which results in 
an effective time of observation of 1173 days and an 
exposure of 2.003-10^'^ m^-s-sr. Simulations show that 
for the above conditions full trigger and reconstruc- 
tion efficiency at KASCADE-Grande is found above 
E « IQi^eV. 

Simulations were employed to study the system- 
atic uncertainties, the performance and reconstruction 
methods in the experiment as well as the influence 
of the cuts. Both the air shower production and de- 
velopment were realized using CORSIKA Q and the 
hadronic generators FLUKA and QGSJET II Q. 
Events with an homogeneous core distribution and 
isotropic arrival direction were produced for a 7 = — 2 
power law spectrum {E — 10^^ — 3 • 10^* eV) and sets 



were generated for a mixed composition of H, He, CO, 
Si and Fe on similar abundances. A weight function 
was finally applied to the simulated data to describe a 
spectrum with 7 = —3, which resembles the observed 
value. 



B. Strategy 

To reconstruct the energy spectrum, three indepen- 
dent techniques have been applied here: one based 
solely on the total charged number of particles 
another one on the muon number [lo| . and the last 
one on a combination of both observables [ll| . 

In the ideal case, when the measurements are ac- 
curate enough, when the reconstruction procedures 
work properly, when the Monte Carlo simulations are 
a faithful description of the EAS and the cosmic ray 
composition is known one expects the same spectrum 
from all three methods. But since these conditions are 
not fulfilled, some differences are expected among the 
final experimental results. At this point, the present 
strategy shows different advantages, it allows (1) to 
carry out different cross-checks of the reconstruction 
procedures, the influence of systematic uncertainties 
and the performance of the detectors, (2) to test the 
sensitivity of Nch and to the elemental composi- 
tion and (3) to study the validity of hadronic inter- 
action models. In the next section, the different re- 
construction methods and results will be described. It 
is worth to say that the reconstruction of the energy 
spectrum has also been performed via the density of 
carged particles at a distance of 500 m [T^ , but the re- 
sults of this study deserves a more detailed discussion. 
They will be presented in an upcoming paper. 
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FIG. 3: The Nch differential spectra derived from 
KASCADE-Grande measurements. 



10" 
10 
10 
10" 
10" 
10" 
10" 
10 
10 
10 





• 0.00° < 


) < 16.71° 




1 * 16.71° <e 


< 23.99° 


O 




< 29.86 


I 


8 a 29.86° <f 


< 35.09° 


^ A 


* . 35.09° <f 


< 40.00° 


i a 














A 






i 






■ 






i 

t 






1 






1 






















F , , , , 1 , 








4.5 



5.5 



6.5 



7 7.5 
log (N ) 



FIG. 4: The A'^^ differential spectra derived from 
KASCADE-Grande measurements. The muon correction 
function was applied to the data. 



C. Reconstruction with a single observable 

The simplest way to reconstruct the energy spec- 
trum is by using one single parameter, as N^h or 

• Data is divided in five zenith angle intervals 
with equal acceptance. In each case, the correspond- 
ing shower size spectrum is extracted (see figures [3] 
andHI and the constant intensity cut (GIG) method 
is applied to correct for attenuation effects in the at- 
mosphere. A reference angle, 9ref, must be chosen for 
this task. In this study, Oref corresponds to the mean 
of the corresponding zenith angle distribution, which 
for Nch is 20° and for 7V^, 22°. The difference comes 
from the particular cuts applied in each case. 

As a next step, from the corrected observable the 
energy is inferred event-by-event by invoking a Monte 
Garlo calibration function, which is composition de- 
pendent. The calibration formula is obtained by fit- 
ting with a power law expression, E = a ■ N^j^^^y the 
Monte Garlo data points for true energy vs shower size 
around 9ref for different mass groups. Then, all data 
from the distinct zenith angle intervals is combined 
according to the energy and a single energy spectrum 
is obtained. 

In a final step, the effect of migration of events in 
the reconstructed energy spectrum is taken into ac- 
count by applying a response matrix, i?y (also ex- 
tracted from simulations). If n1^^ and n*''"'^ are the 
number of events inside the reconstructed and true en- 
ergy intervals, log^Q -E^ and log^Q i^*''"^, respectively, 
then the response matrix is applied in the following 
way: nf "^ = Y.i=i R^s'^T^ ^ where Y.i=i Rij = 1> 
since it represents the logj^g probability distri- 

bution of an event with reconstructed energy inside 
the bin logj^Q-E^. Before unfolding, the response ma- 
trix has to be smoothed to avoid the introduction of 
artificial fluctuations due to the limited statistics in 
Monte Garlo simulations. In figure [5l the unfolded 
spectra are plotted. 



D. Nch - Nf_, method 

This path to reconstruct the energy spectrum ex- 
ploits the information from both the Nch and N^ 
observables The energy of an EAS is derived 

through a Monte Garlo expression, which involves the 
values of the total number of charged particles and 
the muon size in order to reduce the composition de- 
pendence of the energy assignment for the EAS. The 
formula has the form 

logiQ{E/GeV) = [op + {ape ~ ap) ■ k] ■ \ogj^Q{Nch) 
+ [bp + {bpe - bp) ■ k] (1) 

where k — k{Nchi N^) is a mass sensitive parameter 
defined as 

^ ^ logio{Nch/N^) - log,o{Nch/Np^)p 

logi„{Nch/Np)Fe-logio{Nch/Np)p 

in such a way that for protons the average value of k 
is close to zero and approximately one for iron nuclei. 
The constants (a, b) and (c, d) are obtained from fits 
to the scatter plots of E vs Nch and Nch/N^ vs Nch, 
respectively, and performed in the region of maximum 
efficiency. To take into account the influence of the at- 
mospheric attenuation in the EAS the above formulas 
are built for every zenith angle interval and applied 
accordingly. As it was the case for the single param- 
eter reconstruction, the energy spectrum obtained by 
the present method is corrected for the effect of the 
migration of events. The final spectrum is shown in 
figure [H 

E. Flux uncertainties 

Each method has its own systematic uncertainties, 
which were carefully estimated when calculating the 
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FIG. 5: Reconstructed all-particle energy spectrum, 
multiplied by , as obtained from the KASCADE- 
Grande data and the different methods described in the 
present paper. The grey band represents the preliminary 
systematic uncertainty for the spectrum reconstructed 
through the Nch ~ Nf_i method. 



energy spectrum. In the case of the single parame- 
ter reconstruction the sources of systematic uncertain- 
ties considered for this work were the following ones: 
the energy calibration relation, the CIC method, the 
muon correction function, the uncertainty of the spec- 
tral index and the response matrix. For the Nch — 
method the contributions to the systematic errors of 
the spectrum come from the atmospheric attenuation, 
the energy calibration function, uncertainties in the 
primary composition (applying the method to MC 
simulations with different composition assumptions), 
the spectral slope, and the accuracy of the reconstruc- 
tion of shower sizes. The total systematic uncertain- 
ties of the reconstructed spectra at 10^'' eV for the 
different methods are shown in figure [51 



obtained through the Nch approach alone. 

Some general insights about the composition of cos- 
mic rays at high energies can also be extracted from 
these graphs. By departing from the idea that the 
true energy spectrum should lie inside the solutions 
obtained from the N^ and the Nch methods, taking 
into account that the energy spectrum obtained from 
N^ by assuming pure protons (iron nuclei) is higher 
(lower) than the one derived from Nch and looking 
at the position of the spectrum from the Nch — N^^ 
method, then it is noticed that the common solu- 
tion spanned by the different methods favors a rel- 
atively heavy composition at high energies inside the 
QGSJET II/FLUKA framework. 

To investigate the details of the energy spectrum de- 
rived from the Nch — N^ approach, a residual plot was 
constructed by a direct comparison with a flux propor- 
tional to £;~3 0i5 ^ggg figure 0. The spectral index of 
the reference flux was obtained by fitting the middle 
range of the experimental spectrum, i.e. the interval 
E = 10^^'^ — 10"'^'' eV. Two features shows up at figure 
[HI one is a concavity above 10^^ eV and another one is 
a small break at « 10^^ eV. Both structures arc found 
to be statistically significant. For example, a fit with 
a power law spectrum above 10^^ eV gives a spectral 
index of 7 = —3.24 ± 0.08. A preliminary statistical 
analysis with the F-test shows that the significance 
associated with this change of the spectral index is at 
a level of « 99.8%. 

In figure [71 the energy spectrum of cosmic rays ob- 
tained from the KASCADE-Grande measurements is 
compared with the spectra obtained by other cosmic 
ray instruments. In general, a good agreement is ob- 
served at low and high energies with the KASCADE- 
Grande data. 

By now the analysis of the influence of the hadronic 
interaction models is restricted by the statistics of 
the Monte Carlo data sets generated with alternative 



IV. DISCUSSION 

As it was pointed out before, the reconstruction 
of the energy spectrum by using several methods al- 
lows to perform several tests and cross-checks. First, 
by a direct comparison of the resulting all-particle 
energy spectra from the KASCADE-Grande data, a 
good agreement is found inside the respective system- 
atic uncertainties, which means that the experiment 
and its components are well understood, also that 
the reconstruction techniques are working as expected 
and that the hadronic interaction model employed, 
QGSJET II/FLUKA, is intrinsically consistent. 

Another interesting observation from figure [SI it is 
that the muon size at sea level is a very good and com- 
position independent energy estimator compared with 
the total number of charged particles. This conclusion 
is deduced when comparing the difference between the 
solutions derived from the N^^ method against the one 
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FIG. 6: Residual plot for the reconstructed all-particle 
energy spectrum obtained from KASCADE-Grande 
through the Nch — Nf^ method. The systematic error 
band is also shown (dotted lines). 
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FIG. 7: A comparison of the all-particle energy spectrum reconstructed from the KASCADE-Grande data through 
the Nch — Nfj, method with the spectra derived by other experiments. 



hadronic models. However, some preliminary stud- 
ies applying the reconstruction Nch method [13] and 
the two-observables approach have been already per- 
formed with EPOS vl.99 [3]. Those analysis show 
that the differences between the reconstructed spectra 
using EPOS vl.99 and QGSJET II are of the order of 
10—15%. The same preliminary tests done with EPOS 
vl.99 also suggest that the observed stuctures at the 
energy spectrum reconstructed with QGSJET II are 
not an artifact of the hadronic interaction model. 



V. CONCLUSIONS 

The all-particle energy spectrum of cosmic rays was 
reconstructed from the KASCADE-Grande data using 
three different techniques and the hadronic interac- 



tion models QGSJET II and FLUKA. The resulting 
energy spectrum shows statistical significant features 
that can be identified with a concavity and a small 
break at w 10^^ and 10^^ eV, respectively. The na- 
ture of this structures are under discussion. Infor- 
mation about the composition in this energy inter- 
val will provide valuable clues to solve the mistery. 
These composition studies are now right underway at 
KASCADE-Grande. 
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